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Abstract 

Multiparticle production at high energies is described in terms of color strings stretched between the 
projectile and target. As string density increases, overlap among the strings leads to cluster formation. 
At some critical density a macroscopic cluster appears, spanning the entire system. This marks the 
percolation phase transition. Data from p+p, d+Au and Au+Au collisions at 200 GeV has been analyzed 
using the STAR detector to obtain the percolation density parameter, r\. For 200 GeV Au+Au collisions, 
the value of r\ is found to lie above the critical percolation threshold, while for 200 GeV d+Au collisions 
it is below the critical value. This supports the idea of string percolation, which at high enough string 
density is a possible mechanism to explore the hadronic phase transition to a quark-gluon plasma. 
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1 Introduction 

It is postulated that in the collision of two nuclei at 
high-energy, color strings are formed between projec- 
tile and target partons. These color strings decay into 
additional strings via q—q production, and ultimately 
hadronize to produce the observed hadron yields [T]. 
In the collision process, partons from different nucle- 
ons begin to overlap and form clusters in transverse 
space. The color strings are of radius r = 0.20-0.25 
fm [J. 

The fusion of strings to form clusters is an evolution 
of the Dual Parton Model [2J, which utilizes inde- 
pendent strings as particle emitters, to the Parton 
String Model (PSM), which implements interactions 
(fusion) between strings [3]. At some point, a cluster 
will form which spans the entire system. This is re- 
ferred to as the maximal cluster and marks the onset 
of the percolation threshold. An overview of perco- 
lation theory can be found in the following reference 
4 . The quantity T], the percolation density param- 
eter, can be used to describe overall cluster density. 
It can be expressed as 



as [5] 
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with N the number of strings, S the total nuclear over- 
lap area, and tttq the transverse disc area. At some 
critical value of 77 = r/ c , the percolation threshold is 
reached. rj c is referred to as the critical percolation 
density parameter. In two dimensions, for a uniform 
string density, and in the continuum limit, r\ c = 1.175 
0. 

To calculate the percolation parameter, 77, a parame- 
terization of pp events at 200 GeV is used to compute 
the pt distribution 
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where a, po, and n are parameters fit to the data. This 
parameterization can be used for nucleus-nucleus col- 
lisions if one takes into account the percolation of 
strings by [1] 
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In pp collisions at 200 GeV, the quantity 
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1.0 ± 0.1, due to low string overlap probability in 
pp collisions. Once the px distribution for nucleus- 
nucleus collisions is determined, the multiplicity damp- 
ing factor can be defined in the thermodynamic limit 
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which accounts for the overlapping of discs, with 1 — 
e~ v corresponding to the fractional area covered by 
discs. 



2 Data Analysis 

The data utilized in this analysis was acquired by the 
STAR experiment at RHIC. Only data from within 
± 1 unit of pseudorapidity with greater than 10 fit 
points in the main STAR detector, the Time Projec- 
tion Chamber, was used. The events studied were 
minimum bias events, with standard STAR central- 
ity cuts as defined by the trigger detectors and over- 
all event multiplicities, which increase with increasing 
collision centrality. The longitudinal (z) vertex was 
constrained to within ± 30 centimeters of the pri- 
mary event vertex. The distance of closest approach 
(dca) of primary tracks to the event vertex was less 
than 3.0 cm. The transverse momentum of charged 
particles was less than 1.2 GeV. 
The collision systems studied includes: Au+Au at 
200, 62.4, and 19.6 GeV and pp, d+Au, and Cu+Cu 
at 200 GeV. 



3 Results and Discussion 

The percolation density parameter, 77, has been de- 
termined for several collision systems and energies. 
These results have been compared to the predicted 
value of the critical percolation density, rj c . If rj c is 
exceeded it is expected that the percolation thresh- 
old has been reached, indicating the formation of a 
maximal cluster that spans the system under study. 

Figure 1 is a plot of the quantity F(?7) (EqHJ ver- 
sus the percolation density parameter, (77), for central 
collisions. The red line is a representation of F(?7) 
and is drawn to guide the eye. The critical perco- 
lation threshold is indicated as that for overlapping 
strings with uniform density [5|. There is a general 
increase in 77 with increasing system size and energy. 
The predicted value of ?7 C lies just above that of 19.6 
GeV Au+Au collisions and well above that of 200 
GeV d+Au. Both 200 and 62.4 GeV Au+Au lie well 
above the predicted percolation threshold. It is ex- 
pected that the overlap probability for string clus- 
ters increases as the system size and energy increase. 
Therefore, a larger percolation density should be ob- 
tained at 200 GeV Au+Au compared to lower ener- 
gies and other, lighter systems. 
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Figure 1: Multiplicity suppression factor, F(?7) versus the 
percolation density parameter, r\. The line is the function F(rj) 
and is drawn to guide the eye, not as a fit to the points. The 
estimated critical percolation density for 2-D overlapping discs 
in the continuum limit, rj c , is shown. The rj values are calcu- 
lated for the most central collisions. The values for 200 GeV 
d+Au and 19.6 GeV Au+Au lie below rj c , while those for 62.4 
and 200 GeV Au+Au lie above. 



One can also consider the percolation density as a 
function of centrality in Au+Au collisions. The cen- 
trality expressed in terms of the number of partici- 
pating nucleons (N part ) as found from Monte Carlo 
Glauber calculations [7]. More central collisions cor- 
respond to greater values of N z . Figure 2 shows 77 as 
a function of the number of participant nucleons in 
Au+Au collisions at 200 and 62.4 GeV. It is shown 
here that for all centralities, except for the most pe- 
ripheral bin, all 200 GeV Au+Au collisions lie above 
the critical percolation threshold. The value of 77 
increases with increasing collision centrality, an ex- 
pected indication of additional string overlap in more 
central collisions. For 62.4 GeV Au+Au, almost all 
centralities, except for the three most peripheral bins, 
lie above the critical percolation density. Another 
interesting feature is that the value of 77 appears to 
achieve a maximum value of approximately 1.8 in cen- 
tral 62.4 GeV Au+Au collisions, whereas in central 
200 GeV Au+Au collisions 77 continues to increase. 
This could potentially indicate a saturation of the 
cluster overlap area, and hence the percolation pa- 
rameter in 62.4 GeV Au+Au collisions. 

Finally, Fig. 3 presents the fractional area covered 
by clusters as a function of 77. The data is shown for 
all centralities of 62.4 and 200 GeV Au+Au collisions, 
as well as new data with 200 GeV Cu+Cu results. 
The 200 GeV Cu+Cu result covers a region between 
that of mid-central to central 62.4 GeV Au+Au and 
mid-peripheral 200 GeV Au+Au. Because 77 is de- 
pendent on energy and atomic number, this behavior 
is naively expected. Unlike the Au+Au results, the 
Cu+Cu data is not corrected for efficiency and detec- 
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Figure 2: The percolation density parameter, rj, as a function 
of collision centrality (Npart) in 62.4 and 200 GeV Au+Au 
collisions. Larger values of N par t correspond to more central 
collisions. For almost all collision centralities, 200 GeV Au+Au 
exceeds the critical percolation threshold, r\ c . In 62.4 GeV 
Au+Au, all except the three most peripheral bins exceed rj c . 



tor acceptance. It has been predicted by Satz [5] that 
the study of Cu+Cu collisions at RHIC energies pro- 
vides one of the best possibilities to study the onset of 
deconfincment due to the fact that for heavier nuclei 
(such as Pb+Pb) at top RHIC energy, all centralities 
are above the percolation threshold. As shown in Fig. 
3, 200 GeV Cu+Cu and 62.4 GeV Au+Au appear to 
probe the critical percolation threshold. 

4 Summary 

In summary, we have presented the preliminary re- 
sults for the percolation density parameter, 77, at RHIC 
for several collision systems as a function of collision 
centrality. Percolation theory describes the hadronic 
to partonic phase transition as occurring at a critical 
percolation density, -q c . These results show that for 
light systems (d+Au) or heavy systems at low ener- 
gies (Au+Au at 19.6 GeV), the calculated value for 
?7 lies below the predicted critical value. For min- 
bias Au+Au collisions at 62.4 and 200 GeV, 7/ lies 
well above the predicted percolation threshold. As 
a function of centrality in Au+Au collisions at 62.4 
and 200 GeV, almost all centralities except the most 
peripheral for top RHIC energy are above the perco- 
lation threshold, while for the lower 62.4 GeV events, 
several peripheral bins lie just at, or below, r/ c . We 
have also shown preliminary results for Cu+Cu colli- 
sions at top RHIC energy compared to 62.4 and 200 
GeV Au+Au results. The Cu+Cu results span the 
regime between central 62.4 and peripheral 200 GeV 
Au+Au. With increasing 77, the fractional area of 
the system covered by clusters increases. This be- 
havior is expected due to the overall increase in the 
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Figure 3: The fractional area covered by clusters as a function 
of the percolation density parameter, r). Results for 200 and 
62.4 GeV Au+Au, along with 200 GeV Cu+Cu are shown. 
The new, preliminary Cu+Cu results lay in the regime that 
spans the range of r\ covered by peripheral 200 GeV Au+Au 
and central 62.4 GeV Au+Au. 

number of strings as a function of collision energy 
and atomic number of the colliding nuclei. These re- 
sults indicate that 200 GeV Cu+Cu and 62.4 GeV 
Au+Au may provide the ability to study the onset 
of the critical percolation threshold as a function of 
collision centrality. The onset of this collective behav- 
ior in traditional percolation theory is indicative of a 
phase transition, which may provide clues to whether 
a QGP is formed in heavy-ion collisions at RHIC. 
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